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Motivated by the recent emergence of a new class of anisotropic 2D materials, we examine their 
electromagnetic modes and demonstrate that a broad class of the materials can host highly direc¬ 
tional hyperbolic plasmons. Their propagation direction can be manipulated on-the-spot by gate 
doping, enabling hyperbolic beams reflection, refraction and bending. The realization of these nat¬ 
ural 2D hyperbolic media opens up a new avenue in dynamic control of hyperbolic plasmons not 
possible in the 3D version. 


A hyperbolic material (HM) is a highly anisotropic ma¬ 
terial, with real parts of principal components of its per¬ 
mittivity tensor having opposite signs [TJ|2]. The name of 
the material is derived from the topology of the /c-surface 
which is a hyperboloid (unlike an ellipsoid in conventional 
anisotropic materials [3]). Due to this unique topology, 
waves with very large (theoretically infinite) wave vec¬ 
tors can propagate inside HMs. Existence of the high-/c 
modes means very large photonic density of states inside 
the material and thus these materials provide unique ca¬ 
pabilities for controlling properties of light emitters [4], 
creating broadband thermal sources [5], supporting direc¬ 
tional surface plasmons [6l [7], etc. All HMs known so far, 
both natural [8lJl 2 j and artificial mm, are non-tunable 
materials with optical response predefined by the design. 
The hyperbolic regime in metamaterials and metasur¬ 
faces relies entirely on geometrical arrangement of the 
material constituents, and thus is extrinsic in nature. As 
the geometry of a metamaterial is difficult to modify af¬ 
ter the material is created, its optical response at a given 
frequency can not be changed on-the-spot. 

In this letter we argue that the emerging class of 
anisotropic 2D materials [13] is naturally suited to 
serve as HMs capable of supporting tunable hyper¬ 
bolic plasmons. The family of 2D materials is very 
broad (encompassing many elements from the periodic 
table), with tremendous variety in electronic and opti¬ 
cal properties, including metals, semimetals, semicon¬ 
ductors and dielectrics [T4|. Recently, an emerging class 
of 2D materials with anisotropic electronic and opti¬ 
cal properties, such as group V mono- and multi- lay¬ 
ers, most notably black phosphorus (BP)[I5H22], the IT 
phase of transition metal dichalcogenides (TMDs) [23l - 
[25| and trichalcogenides [26] - [28] are garnering attention. 
Even though these materials are semiconductors with a 
bandgap varying from the mid-IR to the optical range, 
injection of free carriers is possible by doping. Doping 
can be controlled by a gate, allowing for maximum car¬ 
rier concentrations up to 10^^ cm~^[29]. 

We offer a novel physical mechanism, inherent in na¬ 
ture, for creating highly-tunable two-dimensional hyper¬ 
bolic materials based on the anisotropic 2D materials in 


the far- and mid- infrared frequency range. This mecha¬ 
nism relies on the interplay between two types of elec¬ 
tron motions defining the electromagnetic response of 
typical semiconductors; motion of electrons (holes) in¬ 
side the conduction (valence) band of the material (in¬ 
traband motion), and electron transitions from valence to 
conduction band of the material (interband transitions). 
In isotropic material predominance of intraband motions 
(interband transitions) renders capacitive (inductive) re¬ 
sponse. We argue, that in novel highly anisotropic 2D 
materials, anisotropic electron motion results in induc¬ 
tive material response along one of the optical axis, and 
capacitive response along the other axis, thus effectively 
rendering a hyperbolic 2D material. As the intraband 
motion of electrons in 2D materials can be tuned on- 
the-spot by applying an electric bias, so does the optical 
response of the resulting 2D hyperbolic material. More¬ 
over, due to its inherent nature, the hyperbolic regime 
can be implemented using a single sheet of 2D material 
without requiring complicated pattering of the substrate. 
We discuss properties of these hyperbolic plasmons and 
examine their existence across the material parameter 
space. 

We begin by considering a minimal model for the con¬ 
ductivity tensor for an anisotropic semiconducting media, 
which accounts for both the intraband electron motions 
(first term on the right hand) and for interband electron 
transitions (second term). 
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j = x,y, where n is the concentration of electrons, rrij is 
the electron’s effective mass along the jth direction, is a 
relaxation time, ujj is the frequency of onset of interband 
transitions for the jth component of conductivity, and Sj 
accounts for the strength of the interband component. It 
may seem somewhat counter intuitive that the x and y 
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components of the conductivity tensor can have different 
frequencies of inter-band transitions. However, this is the 
case in BP multi-layers, where different crystallographic 
directions have different symmetries, and thus the dipole 
moment for electron interband transitions is zero along 
one of the directions [18]. 

We use the Drude model in order to calculate intra¬ 
band conductivity. Interband conductivity is introduced 
phenomenologically by defining absorption through a 
step function The imaginary part of interband 

conductivity then follows from the Kramers-Kronig rela¬ 
tions. In the general case, parameters in Eq. should 
be either measured experimentally or estimated theoret¬ 
ically using, for example, the Kubo formula. Assuming 
that material parameters are isotropic, Eq. can de¬ 
scribe the conductivity of graphene [30|. On the other 
hand, in the anisotropic case this expression can be used 
as an approximation for the conductivity of multi-layer 
black phosphorus [20] . 

We start our consideration with a particular set of 
parameters (see Eig. and then study how the ma¬ 
terial response changes when one or several parameters 
vary. The conductivity of a 2D material is presented in 
Fig. [T] In general, two distinct regimes in the electro¬ 
magnetic response of 2D material can be identified. If 
frequency is sufficiently low (to the left of dash-dotted 
red or blue lines in Eig. Hf), then conductivity is of a 
pure Drude type (Ima > 0) and we refer to the material 
as purely anisotropic (in order to distinguish from hy¬ 
perbolic anisotropic regime). If frequency is sufficiently 
high the contribution from interband electron transitions 
may become dominant, and the imaginary part of con¬ 
ductivity becomes negative. Due to the crystal asym¬ 
metry, Ima^x and Imcr^^ can in principle change sign 
at different frequencies and a frequency interval where 
Imaxx • IiacTyy < 0 exists, which indicates the hyperbolic 
regime (to the right of dash-dotted red or blue lines). 
Such natural hyperbolicity in anisotropic 2D materials 
is unlike hyperbolic metamaterials that have been ex¬ 
tensively studied so farm [6]. Moreover, 2D materials 
offer unique opportunity of switching between the two 
different plasmonic regimes on-the-spot by applying an 
electrical bias (compare two doping cases in Eig. Et)- 

We study properties of plasmons in pure anisotropic 
and hyperbolic regimes by exciting surface waves in a 
2/imx2/im sheet of 2D material by a ^-polarized elec¬ 
tric dipole placed on top of the material (see Eig. 
for details). We choose material parameters to be the 
same as in Eig. ^ Hereinafter we assume that the 
medium above and below the 2D material is vacuum. 
The distribution of plasmon electric field (obtained solv¬ 
ing Maxwells equations numerically using a commer¬ 
cial finite-difference time-domain method (EDTD) from 
Lumerical[3T]i is presented in Eigs. [^-e. Eigs. ^,h 
demonstrate the plasmon field intensity, E, in the purely 
anisotropic regime. In this case the plasmon propagates 
along one of the crystallographic axes of the material {x- 
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FIG. I. Real (a) and imaginary (b) parts of conductivity of 
2D material. Parameters of 2D material are nix = 0.2mo, 
niy = mo, r] = O.OI eV, Wx = 1 eV, Sx = I.7so, Wy = 0.35 eV, 
and Sy = 3.7so. so = jAh, mo is the free-electron mass. 


axis in this case). Eigs. If-e show plasmon propagation 
in the hyperbolic regime. Plasmon energy is channeled 
as very narrow beams, the angle between direction of 
beams propagation and x-axis being a function of elec¬ 
tron concentration n. Thus, the direction of plasmon 
propagation in the hyperbolic regime can be controlled 
on-the-spot by biasing the 2D material. This opens a 
new avenue for solving the problem of plasmon guiding, 
which in conventional materials requires structuring of 
the material geometry |32| or applying external magnetic 
fields j33l [34] . 

In order to understand the behavior of surface waves it 
is instructive to inspect the plasmon dispersion relation 
in 2D materials. Considering only eigenmodes bounded 
to the 2D material plane, 2 : ^ 0), we 

obtain [ 35 ] 
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where p = y Qx ^ ~ ^0 = uyeopo is the vac¬ 

uum wavenumber, and 5o, po are vacuum permittiv¬ 
ity and permeability, respectively. Solution of Eq. (|^, 
q(cc;), defines plasmon wavefront propagation. However, 
in anisotropic materials in the general case directions of 
wavefront and energy propagation do not coincide. In¬ 
stead, it is group velocity Vg = VqCi;(q) that is in the 
direction of plasmon energy propagation [3]. Analytical 
computation of group velocity requires the explicit form 
of co’(q), which is not available in our case. Instead, 
we use /^-surfaces, uj{qx^Qy) = const, for further analy¬ 
sis (see Eig. §)■ Group velocity, and thus direction of 
the plasmon energy flow, has to be orthogonal to the k- 
surface. Assuming that conductivity is purely imaginary 
ajj = ia'Jj (losses are small in our system, see Eig. fTl), 
and that Qx^Qy ^ ^0 5 Eq. © can be further simplifiea 
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FIG. 2. (a)-(e) Spatial distribution (log scale) of electric field, 
|E|, of a surface plasmon excited by y-polarized electric dipole 
(green arrow in panels (a)-(e)) in a 2 x 2 fim sheet of 2D 
material placed in plane z = 0 (dashed black line in panels 
(b,d)). Distribution of electric field is calculated in planes 
z = 10 nm (a,c,e) and x = 0.6 /xm (b,d). (a,b) n — 10^^ 

cm“^, huj — 0.165 eV, (c,d) n — 10^^ cm“^, and hw = 0.3 
eV, (e) n = 3 X 10^^ cm“^, huj = 0.3 eV. (f) k-surfaces, 
oo(qx,qy) = const, for plasmons in panels (a)-(e). 


In the purely anisotropic case {cf'^x^^yy > 0)? 
right side were constant, Eq. Q would be an equation 
for an ellipse in q-space with the axis oriented along q^ 
and qy. The length of the ellipse’s principal axes along 
Qx{Qy) direction is proportional to CFyy{cFxx)‘ Thus, the 
/c-surface always takes quasi-eliptic form elongated along 
the direction of the smallest component of the conduc¬ 
tivity tensor, the degree of elongation being defined by 
the ratio of and ciyy. This is what we see in Fig. 
(blue line), where = 0.55 mS, ciyy = 0.039 mS. Due 
to the strong elongation of the /^-surface along the qy- 
axis, group velocity for predominant range of plasmon 
momentum q points approximately along the qx axis, 
which causes plasmons to carry energy in the direction 
of the X crystallographic axis. If the ratio of and cr'^y 
decreases, the plasmon beam spreads. Finally, in the case 
of (j'l^ = (Tyy the /c-surface turns into circle and there is 
no preferential direction for energy transfer. 

In the hyperbolic case {cf'xx ’ ^yy < 0) ^he /^-surface 
is a hyperbola (see Fig. red and green lines), 

with the hyperbola asymptotes being defined by qy = 

(see Eq. @)- Normals to the asymp¬ 
totes always point along the same direction defined by 

V = ^^\jWyyl<x\ (5) 


which is also the direction of the plasmon beams. As 
plasmons with all momenta carry energy along the same 
direction, the beams are significantly more localized than 
in the purely anisotropic case. Moreover, as can be seen 
from Eq. (§, beams are leaning towards the crystallo¬ 
graphic axis with the largest conductivity. From Fig. 
it is clear that the conductivity of the 2D material, and 
thus the direction of plasmon beam propagation, can be 
changed on-the-spot by doping. Indeed, in the case of 
Fig. 1^ (n = 10^^ cm“^) = 0.3 mS, (jyy = —0.11 

mS, ^ich corresponds to the propagation angle of 31°, 
while in the case presented in Fig. (n = 3 x 10^^ cm“^) 
^xx = 0-076 mS, (Jyy = —0.15 mS, and propagation angle 
is 55°. It should be noted that the propagation angles 
obtained using Eq. (§ agree with the results of numeri¬ 
cal simulations presented in Fig. Finally, we want to 
point out that the beam can experience multiple reflec¬ 
tions from the edges of a 2D material sheet, as can be 
clearly seen from Fig. [^. This leads us to the discussion 
of my optics of the hyperbolic surface plasmons, reminis¬ 
cent of conventional geometrical optics, in the later part 
of the paper. 

It should be noted that there is nothing special about 
the choice of parameters in Figs. |1|2[ Anisotropic and 
hyperbolic plasmons are supported by 2D materials with 
a broad range of parameters. It is quite obvious that as 
long as the material is anisotropic and has free electrons, 
pure anisotropic plasmons can propagate in this material 
when the frequency is sufficiently low. In fact, multi-layer 
BP is an example of such a 2D material that is capable 
of guiding anisotropic plasmons [19]. 

Occurrence of a hyperbolic regime requires an inter¬ 
play between electron intraband motion and interband 
transitions, as it is interband transitions that cause the 
imaginary part of conductivity to change sign. Thus, the 
hyperbolic regime comes into play at frequencies suffi¬ 
ciently close to the frequency of interband transitions. 
In Fig. I^we clarify this notion of closeness by consider¬ 
ing properties of hyperbolic plasmons for a broad range 
of parameters. We assume that rux < ray, which means 
that the x-component of conductivity is predominantly 
Drude-like, while the ^-component of conductivity is the 
one that is mostly affected by interband motions of elec¬ 
trons. 

In Figs. It ,b we inspect the behavior of surface plas¬ 
mons at different frequencies and electron concentrations. 
One can see that the direction of plasmon beams leans 
towards the x-axis at small frequencies and/or high con¬ 
centrations. This is due to large values of a'f^, which 
is predominantly Drude-like (see Eq. compared 

to (Tyy. Moreover, when the frequency is very low and 
electron concentration is very high, intraband motion of 
electrons dominates (jyy as well. This makes the material 
purely anisotropic (thus the cutoff in Figs. [^,b). On the 
other hand, at low concentrations and high frequencies, 
interband motion of electrons is dominant. Conductivity 
\(Jyy\ exceeds \(Jxx\ plasmon beams propagate along 
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FIG. 3. Angle between plasmon beam direction and x-axis, 
for T] = 0.01 eV, rrix = 0.2?tio, rriy — mo, Sx — Sy = 2so, 
hwx — 1 eV, hujy = 0.35 eV, n = 5 x 10^^ cm“^, and huj = 
0.3 eV, unless stated otherwise in the figure. 
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FIG. 4. Plasmon launched by a y-polarized electric dipole in a 
layered medium made by placing sheets of 2D materials with 
different electron concentrations next to each other. Dipole 
position is designated by a green arrow. Parameters of 2D 
materials are the same as in Fig. (a,b) Goncentration of 
electrons in the material above dashed white line is n = 10^^ 
cm“^, while below the white line is n = 3 x 10^^ cm“^. (c) 
Electron concentration varies gradually from n — 10^^ cm“^ 
(upper layer) to n = 2 x 10^^ cm“^ (lower layer). 


directions close to the ^-axis. 

Figs. in ,d demonstrate the behavior of plasmon beams 
at a given frequency and electron concentration, for dif¬ 
ferent values of rrix^ oox and Sy. One can see that any 
change in the set of parameters (decrease of rrix^ increase 
of ujx) that leads to the dominance of intraband motions 
of electrons over the interband ones shifts directions of 
plasmon beams towards the x axis. When ujx becomes 
high, the contribution of interband motions to is neg¬ 
ligible, and further increase of uJx does not change the 
propagation angle, thus the saturation shown in Fig. HI- 
On the other hand, the contribution of interband tran¬ 
sitions to the conductivity increases with an increase in 
Sy, which explains the shift of the plasmon beam prop¬ 
agation direction towards the ^-axis. In general, we can 
see that the hyperbolic regime in 2D materials is robust 
under a very broad range of experimentally achievable 
parameters. 

Due to their outstanding mechanical properties, 2 D 
materials can withstand a significant amount of strain, 
which provides unique capabilities to vary material pa¬ 
rameters in Eq. Q in a wide range using strain 
engineering [36]. Indeed, significant variations in the band 
gap (between 0.38 eV and 2.07 eV) and electron effec¬ 
tive mass (up to 4mo, and relative magnitudes of rrix 
and niy can be interchanged) in BP under the influ¬ 
ence of strain [371439] and/or out-of-plane static electric 
field [ 21 ] were theoretically predicted and experimentally 
demonstrated [40]. Modification of electronic and optical 
properties of ReSe2[24], TiS 3 [ 28 ] under the influence of 


mechanical deformations was also reported. Another po¬ 
tential way of engineering material parameters consists 
in alloying 2D materials or creating heterostructures by 
stacking 2 D materials vertically [dT] |42]. 

The realization of these natural 2D hyperbolic mate¬ 
rials can potentially open up a new avenue in dynamic 
control of hyperbolic surface waves not possible in the 
3D version. In particular, in Fig. we consider prop¬ 
agation of plasmon rays through the interface between 
two sheets of 2D material with different electron concen¬ 
trations placed next to each other. We can clearly see 
reflection and refraction of the plasmons at the interface, 
with reflection angle being equal to the incident angle 
and refraction angle being defined by the conductivity 
of the second medium. Moreover, for a plasmon beam 
incident from the region of higher electron concentration 
(Fig. lit) most of the plasmon’s energy is transferred 
through the interface to the lower doped medium. In the 
opposite case (Fig. most of the plasmon’s energy is 
reflected back to the lower doped medium. Finally, we 
implement an analogue of a graded-index plasmonic ma¬ 
terial by gradually varying the doping across y from 10 ^^ 
cm“^ to 2 X 10^^ cm“^(see Fig. |^). One can see bending 
of the plasmon ray resembling bending of the light ray in 
graded-index medium. 

In conclusion, we demonstrated in a general manner 
for the first time that anisotropic 2D materials can sup¬ 
port propagation of highly directional and tunable hy¬ 
perbolic surface waves. In the purely anisotropic regime 
(Imcr^'^ > 0 , > 0 ) plasmons propagate along the 
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crystallographic axis with larger conductivity, with the 
ratio of the imaginary parts of conductivities defining 
the degree of plasmon localization along the axis. In the 
hyperbolic regime (Imcr^'^ • Ivdcjyy < 0), plasmons propa¬ 
gate as a very localized beam along the direction defined 
by the ratio of the imaginary parts of the conductivi¬ 
ties. Unlike 3D HMs, new 2D natural HMs are highly 
tunable, the frequency range of the hyperbolic regime as 
well as the direction of plasmon beams propagation can 
be modified on-the-spot by gate doping of the 2D ma¬ 
terial. Through doping control reflection, refraction and 
bending of plasmon beams can be achieved. By elucidat¬ 
ing the existence of an interesting hyperbolic regime due 
to the interplay of intra- and inter-band electron motion, 
we hope this will inspire and guide their experimental 
search and verification. 


* tlow@umn.edu 
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